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ABSTRACT 

We  extend  the  results  of  our  previous  work  (Part  I)  concerning  steady 
ideal  flow  with  vorticity  in  a  bounded  domain  to  the  situation  where  the 
domain  is  unbounded  and  the  flow  is  uniform  at  infinity.  Prototypical 
examples  include  vortex  pairs  and  Poppl  vortex  wakes  behind  a  cylinder  (in 
uniform  translation).  The  existence  of  solutions  and  the  asymptotic  behavior 
of  these  solutions  in  a  certain  singular  limit  are  established  as  in  Part  I 
using  a  direct  variational  method.  The  variational  principle  needed  here  is 
rather  nonstandard,  and  so  a  detailed  discussion  of  its  formulation  is 
given.  Special  difficulties  arise  for  flows  in  an  unbounded  domain  due  to  a 
lack  of  compactness  (in  the  appropriate  function  space);  consequently,  we  find 
that  there  is  nonexistence  of  solutions  in  sor  is  cases. 
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SIGNIFICANCE  AND  EXPLANATION 


The  purpose  of  this  sequel  (to  Part  I)  is  to  apply  the  results  detailed 
in  the  previous  work  for  a  simple  model  case  to  a  class  of  more  complicated 
problems  arising  in  the  dynamics  of  a  two-dimensional,  ideal  fluid.  The 
theory  of  steady  flow  with  vorticity  developed  in  Part  I  is  now  extended  to 
include  certain  flow  geometries  commonly  occurring  in  practice  -  that  is, 
steady  ideal  flows  representing  vortex  pairs  or  (symmetric)  vortical  wakes 
behind  a  symmetric  obstacle  in  a  uniform  stream.  As  before,  special  emphasis 
is  given  to  flows  with  concentrated  vorticity.  A  detailed  discussion  of  the 
physical  principles  underlying  the  formulation  of  the  rather  novel  variational 
approach  to  these  problems  is  presented.  Much  of  the  analysis  concerning  the 
qualitative  description  of  solutions  parallels  that  of  Part  I,  although  some 
important  differences  must  be  dealt  with  here. 
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ON  STEADY  VORTEX  FLOW  IN  TWO  DIMENSIONS,  II 


Bruce  Turkington* 

Department  of  Mathematics 
Northwestern  University 
and 

Mathematics  Research  Center 
University  of  Wisconsin-Madison 

In  this  sequel  to  Part  I  we  continue  the  study  of 
steady  ideal  fluid  flows  in  two  dimensions  which 
possess  vorticity.  Specifically,  we  consider  steady 
vortex  flows  imposed  upon  an  underlying  (non-trivial) 
irrotational  flow.  Classical  examples  of  such  flows 
include,  for  instance,  the  uniformly  translating  vortex 
pair  and  the  Fftppl  vortex  wake  behind  a  cylinder  in  a 
uniform  stream;  it  is  the  appropriate  generalizations 
of  flows  such  as  these  that  we  examine  in  the  present 
paper.  As  in  Part  I  we  characterize  solutions  of  the 
Euler  fluid  dynamical  equations  variationally  as 
extremals  for  an  appropriate  energy  functional  of  the 
vorticity.  This  approach  permits  us  to  develop  here 
the  natural  extensions  of  both  the  existence  and  the 
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asymptotic  results  given  in  Part  I.  In  the  present 
context,  however,  the  variational  problem  itself  is 
less  familiar  and  therefore  its  formulation  is 
discussed  in  detail.  Particular  emphasis  is  placed  on 
unbounded  fluid  domains  since  special  difficulties  (not 
all  of  which  are  resolved  here)  arise  in  this  case;  for 
example,  complementary  existence  and  nonexistence 
theorems  are  proved  for  the  vortex  pair  problem. 

§1.  VARIATIONAL  PROBLEM 

Let  D  c_  R  be  the  closure  of  a  bounded,  simply- 
connected  subdomain  of  R+  =  (x=(x^,x2)  6  R^:x2  >  0} 
for  which  0  e  3D.  We  shall  assume  by  suitably  scaling 
the  x  variables  that  5  c  (|x|  <  1).  We  shall  study 
the  flow  of  an  ideal  fluid  with  unit  density  in  the 
(fluid)  domain  defined  by 

(1.1)  D  =  (x  e  R|  :  IxjJ  <  ax,  x2  <  a2)\D 

for  some  (fixed)  truncation  parameters  1  <  a1#a2  <  •. 
We  shall  assume  that  3D  n  R+  is  smooth,  so  that  3D 
is  piecewise  smooth  and  has  no  re-entrant  corners.  We 
refer  to  D  as  the  obstacle.  The  flow  geometry  we 
have  in  mind  is  depicted  in  the  figure. 

The  role  of  the  truncation  parameters 
a  **  (a1#a2)  is  principally  a  technical  one.  It  is  of 
considerable  interest  (and  in  some  ways  more  natural) 
to  study  flows  in  the  unbounded  domain  D  ®  R+\D; 


however,  certain  aspects  of  the  subsequent  analysis  are 
complicated  in  this  case.  Therefore,  we  choose  to 
restrict  our  attention  to  truncated  domains  of  the  form 
(1.1)  while  developing  the  general  results  (§2-4).  We 
then  give  more  specialized  results  for  flows  in 
unbounded  domains  of  the  above  type;  in  §5  flows  past 
an  obstacle  D  *  it  are  treated,  while  in  §6  the  vortex 
pair  problem  for  which  D  =  R+  is  studied.  We  remark 
that  since  3D  n  (x2  =  0}  shall  always  be  a  stream¬ 
line,  all  of  these  flows  may  be  thought  of  as  extended 
(by  reflection)  to  the  symmetric  domain  D*  formed  by 
the  union  of  D  with  its  reflection  about  the  x^  axis. 

The  discussion  in  much  of  the  sequel  parallels 
that  already  given  in  Part  I  (namely  [4]),  and  the 
notation  is  kept  consistent  with  that  of  Part  I.  We 
shall  proceed  to  formulate  the  variational  problem  to 
be  studied  in  the  present  Part  II,  being  brief  at  those 
points  where  the  relevant  details  have  been  developed 
in  Part  I. 


The  (fluid)  velocity  field  is  given  by 

(1.2)  u  =(u1(x),u2(x))  =  JVi|»(x),  x  e  D 

for  some  streamfunction  i>;  recall  J(a^,a2)  =  (a2,-a1). 
The  vorticity  is  then 

(1.3)  u»  =  cd(x)  =  -Aiji(x) ,  xeD. 

We  seek  flows  satisfying  the  following  boundary 
conditions : 


♦  = 

-u 

on 

3D 

*  = 

-Qx2  -  u 

on 

3D\3D 

for  constants  -*  <  u  <  08  and  0  <  Q  <  •.  This  means 
in  physical  terms  that  we  impose  a  flux  across  the 
boundary  segments  3D  n  =  ±  a.)  with  its  rate 
prescribed  by  Q  (the  stream  speed)  while  we  require 
that  the  flow  be  tangential  on  the  remainder  of  3D 
(as  ♦  is  constant  there).  The  value  of  ♦  -  defined 

only  up  to  an  additive  constant  -  is  normalized  by  the 
parameter  u.  Clearly,  these  boundary  conditions  are 
chosen  so  that  a  flow  in  the  truncated  domain  D 
approximates  an  unbounded  flow  past  the  obstacle  D 
with  velocity  (-Q,0)  at  infinity. 

The  Green  function  for  -A  in  D  is  denoted  by 

(1.5)  g(x, x' )  =  log  I x  -  x'j"1  -  h(x,x'),  x,x'  6  D  , 
and  the  Green  operator  by 

(1.6)  Gw(x)  =  /  gfX'X'iufx'idx'  • 

D 


Also,  we  write 


x  e  d  . 


(1.7)  H(x)  »  j  h(x, x) , 

Let  n  *  n(x),  x  e  d,  denote  the  (unique)  solution  of 

An  ■  0  in  D 

(1.8)  n  *  0  on  3D 

n  “  x2  on  3D\3D  ; 

then  n  is  the  streamfunction  for  a  normalized  irrota- 
tional  flow  in  D.  Now  the  streamfunction  ♦  can  be 
expressed  in  terms  of  the  triple  o>,Q,y,  since  accord¬ 
ing  to  (1.3)  and  (1.4)  there  must  hold 

(1.9)  t|»»Ga)-Qn-ji. 

As  explained  in  Part  I#  the  steady  Euler  fluid 
dynamical  equations  can  be  expressed  in  weak  form  as 

(1.10)  /  o>3(*,*)dx  -  0  for  all  *  e  c"(D)  y 

D  u 

recall  3(^,4)  If  Q  is  prescribed  then 

this  is  a  condition  for  u»  alone.  In  fact,  this  is 

the  condition  that  u>  be  a  constrained  extremal  for  a 

certain  (energy)  functional  Eq.  Let  the  vortex  energy 

be  defined  by  the  functional 

(1.11)  E(o> )  “  i  /  u>(x)Gw(x)dx 

Z  D 

=  4  /  /  g(x,  x'  )u>(x)t«)(x,  )dxdx*  , 

D  D 

and  let  the  vortex  impulse  be  defined  by  the  functional 

P(w)  =  /  n(x)u(x)dx  . 

D 


(1.12) 


For  any  given  0  <  Q  <  •  we  consider  the  functional 


(1.13)  Eq(o»)  =  E(«)  -  QP(«)  . 

Then,  as  shown  in  Part  I,  the  dynamical  condition 
(1.10)  is  equivalent  to 


d 

dt 


=  0  ? 


the  variations  of  w  are  defined  by 

w(t)(x)  =  «(£“1(x)),  | t |  <  t,  where  x  =  £t(y)  is 
the  measure  preserving  dif feomorphism  of  D  defined  by 
solving  the  equations 

||  =  JV*(x),  x(0)  =  y  . 

Thus  we  base  our  variational  approach  to  the  problem  on 
the  fact  that  any  extremal  <■»  of  Eq  over  a  class  of 
admissible  functions  containing  all  the  variations 
yields  a  solution  of  (1.10)  -  that  is,  a 
dynamically  possible  flow. 

The  physical  interpretations  of  E  and  P  are 
discussed  in  detail  in  the  appendix.  Only  in  the 
special  case  of  vortex  pairs  in  which  D  *  R+  and 
n(x)  =  X2  do  these  formulas  seem  to  appear  in  the 
literature?  for  example,  see  Batchelor  [1].  Of  course, 
the  corresponding  concepts  of  energy  and  impulse  are 
well  developed  in  the  standard  literature  for 
irrotational  flows. 


We  let  the  class  of  admissible  functions,  K^(D), 
be  those  0  <  u>  e  L*(D)  satisfying  the  constraints 

(1.14)  /  w(x)dx  -  1 

D 

(1.15)  ess  sup  u(x)  <  X  . 

xSD 

The  circulation  in  (1.14)  is  normalized  to  be  1  by 
scaling  u,  and  the  vortex  strength  parameter,  X,  is 
required  to  satisfy  X  meas.  D  >  1.  In  §2  we  give  the 
basic  existence  result  which  states  that  there  exists  a 
maximizer  u  s  uq  x  for  eq  over  the  class  K^(D). 
Furthermore,  the  maximizer  w  has  the  special  form: 

(1.16)  (o  -  Xlfl,  a  -  {x  6  D  i  ( x )  >  0) 

where  <|»  is  defined  by  (1.9)  and  p  is  a  constant 
(uniquely  determined  by  u).  There  is  no  corresponding 
uniqueness  result;  indeed,  nonuniqueness  is  present  in 
some  geometries. 

In  §3-4  we  turn  our  attention  to  the  limiting 
behavior  of  solutions  as  X  ♦  •.  We  show  that 

(1.17)  wQ  A(x)  ♦  5(x  -  Xq)  as  X  +  • 

in  the  sense  of  distributions;  6(z)  is  the  unit 
(Dirac)  delta  measure  at  z  *  0.  The  limit  point  Xq 
is  characterized  by  the  property 

(1.18)  hq<xq)  *  min  Hq(x),  Hq(x)  ■  H(x)  +  Qn(x)  . 

xeD 

The  precise  nature  of  the  limiting  form  of  solutions  is 


expressed  in  terms  of  their  scaled  versions: 


-1. 


u-19>  5Q,x (y>  -  *  “Q,x(xQ,x  +  £y) 


=  i. 


Q»  1 


(Aq  ^  C_  R^  is  open 


and  bounded) 
.2 


where  X0  .  =  f  x  u  .(x) dx  and  Xire^  =  1.  We  prove 

W  t  A  U  #  a 


that  as  X  ♦  • 

^0  2 

L  i  ♦  IR  /rn  weakly  star  L  (R  ) 

(1.20)  •  Q'  l1  1 

3Aq  ^  ♦  dB^fO)  C1  sense  as  curves  . 

The  methods  used  in  §2-4  are  exactly  those  already 
developed  in  the  corresponding  sections  of  Part  1; 
consequently,  these  sections  of  the  present  Part  II  are 
brief. 


New  difficulties  fundamentally  different  from 
those  in  Part  I  arise  in  §5-6  when  unbounded  domains 
are  considered.  Since  then  the  functional  EQ  is  no 

GO 

longer  weakly  (star  L  )  continuous  on  K^(D)  we  are 
forced  to  find  solutions  on  truncated  domains  and  to 
supply  a  priori  estimates  for  their  support.  Only 
under  the  restriction  that  X  be  large  enough  are  we 
able  to  carry  out  such  a  procedure  for  the  general  case 
of  a  domain  D  =  R+\D;  the  previously  established 
asymptotic  properties  of  solutions  are  used  crucially 
in  the  a  priori  estimates  needed  for  the  existence 
proof.  This  analysis  constitutes  §5.  These  steady 
vortex  flows  are  the  appropriate  generalizations  of  the 
classical  example  due  to  Fbppl  of  a  point  vortex  pair 


stationary  with  respect  to  a  moving  cylinder  (in  the 
present  notation  D  -  (|x|  <1,  x2  >  0}  and  X  =*  «); 
the  Fttppl  vortices  are  referred  to  by  Larnb  [3]  p.  223. 
The  families  of  such  flows  (as  constructed  in  §5)  are 
of  some  real  interest  as  they  provide  simple  models  for 
a  familiar  physical  experience*  the  steady  (symmetric) 
vortical  wake  behind  a  symmetric  obstacle  in  a  uniform 
stream. 

The  more  tractable  problem  of  (uniformly 
translating)  vortex  pairs  in  the  special  domain 
D0  ■  R+  is  taken  up  in  §6.  In  this  case  the  vorticity 
possesses  an  additional  symmetry*  a»(-xlfx2)  *  wtx^Xj) 
and  o)(x1#x2)  is  monotonic  in  Xj^  for  xj^  >  0. 

Various  manifestations  of  this  permit  a  more  complete 
analysis  of  solutions  than  is  possible  in  §5.  A 
complete  existence  theory  can  be  given  but  it  requires 
that  the  class  K^(DQ)  be  slightly  widened.  In 
particular,  the  nonexistence  of  solutions  (of  the  form 
(1.16)  for  a  prescribed  Q  >  0)  in  K^(DQ)  is 
demonstrated  for  small  enough  X;  in  terms  of  the  flow 
this  means  that  the  vortex  core  cannot  be  too  broadly 
spread.  Finally,  an  alternate  variational  approach  is 
sketched  which  closely  resembles  the  approach  used  for 
vortex  rings  (the  axisymmetric  analogue)  by  Friedman 
and  Turkington  [2], 


§2.  EXISTENCE 


Let  D  be  as  in  (1.1)  with  a1# a2  fixed.  Let 
K^(D),  the  class  of  admissible  functions  in  D,  be 
defined  by 

(2.1)  K.(D)  =  {w  6  L*(D)  s  /  <o(x )dx  ■  1, 

A  D 

0  <  «(x)  <  X,  a.e.  x  6  D)  . 

Let  the  (energy)  functional  Eq  be  defined  on  K^(D)  by 

(2.2)  E_(w)  =  i  /  /  g(x#x')<i>(x)«(x') dxdx ' 

u  *  D  D 

-Q  /  n(x)a>(x)dx 
D 

for  any  given  0  <  Q  <  •;  recall  the  notations 
established  in  (1. 5-1.8).  We  assume  throughout  the 
sequel  that 

(2.3)  X  >  (meas.  D)”1  . 

The  following  existence  theorem  provides  an 
absolute  maximizer  for  Eq  in  the  class  K^(D). 

Theorem  2.1.  There  exists  “  “  wq  x  e  such  that 

(2.4)  E  (w)  =  max  E  (u)  . 

Q  u6Kx(D)  Q 

Proof.  The  proof  is  almost  identical  to  the  proof  of 
Theorem  2.1  in  Part  I,  except  that  E  there  is 
replaced  by  Eq  here.  We  leave  the  easy  modifications 
to  the  reader. 
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The  computations  sketched  in  §1  (their  more 
detailed  counterparts  are  given  in  Part  I)  yield  the 
following  variational  conditions. 

Corollary  2.2.  Whenever  «•»  e  K^(D)  satisfies  (2.4) 
then  there  holds 

(2.5)  /  u3(Gw  -  Qn, $ )dx  -  0  for  all  ♦  €  c"(D)  • 

D 

The  above  functional  dependence  expressed  in  weak 
form  can  be  given  explicitly. 

Corollary  2.3.  Whenever  «  e  K^(D)  satisfies  (2.4) 
then  there  exists  a  (uniquely  determined)  constant  u 
such  that 

(2.6)  u  *  Xljj  a.e.  in  D, 

ft  -  {x  e  D  s  Gai(x)  -  Qn(x)  >  v }  i 
Ifl  denotes  the  characteristic  function  of  ft. 

Proof.  Again,  the  proof  is  almost  identical  to  the 
proof  of  Corollary  2.3  in  Part  I.  Using  the  former 
methods  we  now  find  that 

u  **  ess  supCGw(x)  -  Qn(x)3  -  ess  inf[Gw(x)  -  Qn(x)J  . 
w(x)<A  w(x)>0 

As  before  we  then  show  that  u  -  0  a.e.  in 
{Gu>  -  Qn  *  p).  We  note  that  in  contrast  to  the  case 
in  Part  I  v  is  not  necessarily  positive. 
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§3.  ASYMPTOTIC  ESTIMATE 


We  let  <u  =  o>q  ^  be  any  maximizer  as  in  Theorem 
2.1,  and  let  ^  be  given  by  (1.9). 


Lemma  3.1.  There  holds 


(3.1) 

V“> 

*  It  1o9  t  -  ci 

(Axe2  =  1) 

where 

is  a 

(positive)  constant 

depending  on  Q. 

Proof. 

We  let 

«  6  K^(D)  be 

(3.2) 

“ =  XIMV 

where  is  chosen  such  that 

H  (X  )  =  min  Hq ( x ) ,  HQ(x)  =  H(x)  +  Qn(x)  . 
u  u  x€D 

Then  a  straightforward  calculation  yields 

V“)  »  eq(“>  >  57  h  '  VV  +  o(1> 

as  A  -*•  «,  and  so  (3.1)  follows. 

The  basic  asymptotic  estimate  for  large  A  is 
proved  by  adapting  the  methods  of  Part  I. 

Theorem  3.2.  There  is  a  constant  R  >  1  depending 
on  Q  (independent  of  A)  such  that 

2 

(3.3)  diam(supp  u> )  <  Re  (Aire  =1). 

Proof.  We  define  as  in  Part  I 

(3.4)  T(uj)=i  /  'J'wdx  ; 

D 
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then  we  have  clearly  an  identity  and  an  estimate: 

(3.5)  E(w)  =  T(u>)  +  j  QP  ( t*» )  +  , 

(3.6)  y  >  2Eq(w)  -  2T(o»)  . 

We  claim  that  T(w)  <  C2  for  a  constant  C2 

depending  on  Q  (independent  of  X).  First  we  notice 

that  y  >  -0/2(3^^  -  1 )  X ;  this  follows  since  if  y  <  0 

then  (x  :  l<lx1l<a1,  Qx2 < I y I >  c  {xeD  :  Qn(x)  <  |y|) 

c  supp  w,  and  hence  2(a]^  -  l)|y|/Q  <  1/X.  Therefore, 

taking  X  sufficiently  large  we  may  assume  that 

y  >  -  i.  Now  we  estimate 

2T(o>)  <  /  (♦-  l)+udx  +  1  ; 

D 

but  since  ♦  <  -y  <  i  on  3D  we  may  apply  the 

reasoning  of  Lemma  3.2  in  Part  I  to  conclude 

/  |V(^»  -  l)  +  |2dx  **/(♦-  1  )+wdx  <  c  , 

D  D 

and  so  the  claimed  bound  for  T(o>)  is  proved. 

Returning  now  to  (3.6)  we  have  for  arbitrary 
x  e  supp  to 

Gu(x)  >  Gw(x)  -Qn(x)  >  y  log  j  -  C3  , 

by  virtue  of  Lemma  3.1.  Hence  we  conclude 

-2irC3  <  /  log  ^tj  <*»(x')dx'  . 

From  this  point  on  the  proof  is  completed  exactly  as  in 
the  proof  of  Theorem  3.3  in  Part  I. 


Remark.  As  is  evident  from  the  above  proof  the 
constant  R  depends  on  the  diameter  of  D,  and  so  the 
resulting  estimate  cannot  be  applied  directly  to 
solutions  defined  in  unbounded  domains. 

§4.  LIMITING  BEHAVIOR 

The  precise  limiting  form  of  w  =  x  as  A  ♦  • 
is  discussed  in  this  section.  By  extracting  a  sequence 
X  «  Aj  -*•  •  if  necessary  we  suppose  that 

(4.1)  XQ  x  =  /  xwQ  X(x)dx  ♦  xj  e  D  as  X  ■  A^  ♦  •  . 

Theorem  4.1.  Any  Xq  as  in  (4.1)  satisfies 

(4.2)  Hq(x£)  -  min  Hq(x),  Hq(x)  =  H(x)  +  Qn(x)  . 

xeo 

Proof.  We  follow  the  method  of  proof  of  Theorem  4.3  in 

A 

Part  I.  Let  u  e  KX(D)  be  defined  by  (3.2).  Then 
(recalling  Lemma  4.2  in  Part  I)  we  have 

E  (<■>)  <  Eq(o>)  <  i-  f  f  log  |  x-x '  |""1co(x)<i>(x'  )dxdx ' 
u  u  D  D 

-iff  h(x,x*  )u>(x)u>(x'  )dxdx'  -  Q  /  n(x)«(x)dx. 

C  D  D  D 

Taking  X  ♦  •  and  using  Theorem  3.2  we  find,  as 
required, 

HQ(X^)  <  HQ(XQ)  -  min  Hq ( x )  . 
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We  are  now  able  to  assert  that  as  X  ■  A^  ♦  • 

(4.3)  u)q  ^(x)  *  “  xq)  t^ie  sense  °f  distributions 

where  5(x  -  Xq)  is  the  unit  (Dirac)  delta  measure  at 

x  *  given  by  (4.2).  Moreover,  the  statements  of 

Theorem  4.4  and  Theorem  4.5  in  Part  I  hold  without  any 

change  for  the  scaled  versions  of  ^  and  ; 

respectively;  we  will  not  restate  them  here.  We  note 

*  * 

that  the  limiting  functions  C  and  v  are,  of 
course,  independent  of  Q. 

§5.  UNBOUNDED  DOMAINS 

In  this  section  we  extend  the  theory  previously 
developed  in  domains  of  the  form  (1.1)  to  unbounded 
domains 

(5.1)  D  «  rJ\d  (D  *  ^) 

where  D  is  as  before;  that  is,  we  study  flows  in  the 
(full)  exterior  of  a  symmetric  obstacle.  We  now  write 
D0  =  R+,  D^  =  R+\{|x|  <  1),  and  we  assume  (as  before) 
that  Di  D  _c  Dq.  Also,  we  denote  by  g(x,x'), 
g0(x,x'),  g1(x,x')  the  Green  functions  corresponding 
to  D,  Dq,  D^  respectively.  The  streamfunction  n 
now  represents  irrotational  flow  past  D  with  velocity 
(1,0)  at  infinity,  and  is  defined  by 
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An  =  0  in  D 
n  =  0  on  3D 

(5.2)  n  o 

n  =  x2  +  o( I x I  x),  Vn  =  (0,1)  +  o( I x I  2) 

as  | x |  -*■  «  . 

The  following  lemma  provides  an  a  priori  bound  for 
the  limiting  location  of  a  solution  as  A  -*•  «*. 

Lemma  5.1.  Let  x*  e  D  be  an  arbitrary  critical  point 
of  HQ  =  H  +  Qn,  0  <  Q  <  «.  Then 

(5.3)  |xj|  <  AQ"1 

(5.4)  xj  <  AQ-1 

for  a  constant  A  depending  only  on  D. 

Proof.  The  proof  of  these  estimates  requires  some 
technical  asymptotic  properties  of  the  functions  H 
and  n,  so  we  present  the  demonstrations  of  these 
properties  first. 


Property  of  H;  Let  H  be  defined  by 


(5.5)  H(x)  «  ^  log  ~~  +  H(x), 


x  6  D  . 


We  claim  that 

15^^  I  57u  /  \  I  e.  r*  I  «l“4 


I 7H ( x ) I  <  Cx2lxJ 


( I x I  >  2)  . 


To  prove  this  we  observe  that  VH(y)  =  7xh(y,y),  y  6  D, 


if  we  define 


h(x,x')  =  h(x,x')  -  ■=—  log  |  x  -  x' 


-1 


where  we  write  x'  ■  (x{,-X2).  For  any  fixed  x'  e  D, 

A  A 

h(x,x')  is  harmonic  in  x  e  D  and  h(x,x')  *  0  when 

A 

x2  *  0.  Thus,  h(x,x')  can  be  extended  by  reflection 
to  be  a  harmonic  function  in  x  e  D*  ■  R2 \(D  u  D  ) 
where  D  is  the  reflection  of  D  about  the  x^-axis; 

A 

h(x,x')  is  then  an  odd  function  of  x2.  It  follows 
that  the  gradient  is  estimated  by 

(5.7)  |Vxh(y,y)|  <  Clyl"1  sup  |h(x,y) I  ( I y I  >  2). 

I  x-yl  <5-lyl 

In  order  to  bound  the  right  hand  side  of  this 

A 

inequality  we  introduce  h^(x,x')  defined  to  be  the 
corresponding  function  for  the  domain 
Di  “  R 2 \ ( I x |  <  1);  it  is  computed  explicitly  to  be 
hj.  ( x,  x 1 )  -  log(l  +  4x2x£|x*  I  “2 1  x  -  x'*f”2)  , 

£  a 

where  we  write  x'  »  |x'|  x*.  Now  since 

g1(x,x>)  <  g(x, x* )  <  g0(x,x')  for  x,x'  e  D1#  we  find 

A  A 

that  0  <  h(x,x')  <  h1(x,x')  for  x,x'  e  Dj^. 

Therefore,  the  desired  estimate  (5.6)  follows  from 

(5.7)  as  we  verify  that 

A  A 

sup  hL(x,y)  <  Cy2 I y I . 

I x-y | <jlyl 

A 

Properties  of  n:  Let  n  be  defined  by 

A 

(5.8)  n(x)  »  x2  +  n(x),  x  e  D  . 

A  A 

We  intend  to  determine  expansions  for  n  and  Vn  for 
large  x?  we  proceed  as  usual  from  the  Green 
representation  formula 


A  jg 

(5.9)  n(x)  =  /_  [g0(x,x')  yyr(  x*  )-n  (x ' )  yyr(x,  x' )  ]ds ' , 

valid  according  to  (5.2).  Expanding  g0(x,x'),  I x I  >  1, 
in  a  Taylor  series  (up  to  second  order)  in  x*  about 
x'  *  0  we  obtain 

(5.10)  n(x)  =  2C2x2lx|"2  +  4C12XiX2 I x|"4  +  0(|x|-3) 
with 

C2  “  ^  3D  ^*2  ^  (X<)  "  nU,)  ^d8‘ 

and  C12  given  similarly  (x£x2  replacing  x2  in  the 
integral);  the  other  corresponding  coefficients 
cl,cll,c22  are  a11  *ero*  Analogously,  first 
differentiating  (5.9)  and  then  expanding,  we  get 

(5.11)  Vn(x)  =  V(2C2x2Ix|“2  +  4CJ2XJX2 | x|”4)+0( |x|“4)  . 

We  claim  that  C2  <  0  (provided  D  *  ^).  Indeed, 
using  the  fact  that  n  =  -x2  on  3D,  we  find 

*  3X 

(5.12)  -2xC2  =  J_  [-X2  yj  +  n  j—lds 

3  D 

»  an  3 x. 

=  /g  Tv  "  x2  yHds 
=  /  |Vn|2dx  +  /  |Vx,|2dx 

D  5 

=  M  +  meas.  D 
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where  M  *  /  I V ( n  -  x2)|2  is  the  (so-called)  induced 
~  D 

mass  of  D. 

The  particular  consequence  of  these  derivations 
needed  in  the  subsequent  proof  is  the  expans ion * 

(5.13)  nx  (x)  -  -4C2x1lxr440(|*r4)  as  I xx  I  ♦  •». 

2  1 

Proof  of  (5.4).  By  the  above  we  have  as  x2  ^  * 

VH(x)  -  (0,  +  7H(x),  |  VH  ( x )  |  -  0(|x|“3)  , 

vn(x)  -  (0,1)  +  Vn(x),  | Vn (x) I  -  o(|x|"2)  . 

Together  these  give 

(5.14)  VH(x)-Vn(x)  -  +  0(|x|"3)  as  x2  ♦  •  . 

2 

But  at  a  critical  point  x  *  x*  there  holds 

VH(x*)*Vn(x*)  »  -Q|Vn(x*)|2  ■  -Q  +  o( | x | ”4)  7 
clearly,  this  contradicts  (5.14)  if  x2Q  >  A  for  a 
sufficiently  large  constant  A. 


Proof  of  (5.3).  For  this  part  we  calculate 


7H(x)»jVn(x)  »  ■  nx  (x)  +  7H(x) »J7n(x)  . 

2  1 


As  IxjJ  *  «  for  x2  <  AQ_1  we  get,  using  (5.13), 
(5.6)  and  (5.2), 

(5.15)  7H(x)‘J7n(x)  =  -  i  C2Xllxr4  +  0(Q_1|xr4)  . 

But  at  a  critical  point  x  *  x*  there  holds 

7h(x  )*J7n(x  )  *  -Q7n(x*)»J7n(x  )  «  0  ; 
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this  contradicts  (5.15)  if  Ix^lQ  >  A  for  a 
sufficiently  large  constant  A,  since  C2  *  0. 

In  the  next  theorem  we  find  that  maximizers  over  a 
truncated  class  yield  solutions  in  the  (entire)  domain 
O  provided  that  X  is  large.  Let 
Da  =*  D  n  {|xjJ  <  a1#  X2  <  a2)  and  let  the  class 
K®(D)  consist  of  all  those  a>  6  K^(D)  such  that 
<0  =  0  a.e.  in  D\Da. 

Theorem  5.2.  Let  a  =  (a1#a2)  be  fixed  so  that 
al'a2  >  ^  max{l,AQ  with  A  as  in  Lemma  5.1.  Then 
any  (absolute)  maximizer  “  “  “q  x  f°r  EQ  over 
K*(D)  satisfies 

(5.16)  supp  wcdxjJ  <  -  1,  0  <  x2  <  aij  -  1) 

provided  X  >  A 

where  A  is  a  sufficiently  large  constant  depending  on 
a,Q,D. 


Proof .  The  existence  of  a  maximizer  w  follows  by  the 
proof  of  Theorem  2.1.  The  basic  asymptotic  estimate 
(3.3)  with  a  constant  R  depending  on  a,Q,D  holds 
for  any  such  maximizer  u>;  the  different  method  of 
truncation  used  here  in  §5  does  not  affect  the  validity 
of  the  proof  given  in  §3.  Taking  A  large  enough  so 
that  X  >  A  implies  Re^  <  and  lXQ  ^  -  XqI  < 


-20- 


(recalling  (4.1)  which  is  valid  in  the  present  case 
also)  we  find  that 

supp  o>  c  Bj^X^) 

for  some  minimum  point  Xq  of  Hq.  Now,  the  required 
containment  (5.16)  follows  immediately  from  Lemma  5.1. 

The  purpose  of  (5.16)  is,  of  course,  to  ensure 
that  the  free  boundary  belonging  to  a  maximizer  w 
does  not  touch  the  truncating  boundary  where  x^  *  ±a^ 
or  x2  “  a2*  It  is  possible  then  to  proceed  as  in  $3 
to  the  conclusion  that  for  »  given  by  Theorem  5.2 
there  holds 

0  *  /  a)3  (Gw  -  Qn,+  )dx  for  all  t  e  c!!(D)  * 

D 

and  (for  a  uniquely  determined  constant  u) 
o)  *  Xlfl  a.e.  in  D,  ft  ■  {x  e  D  :  Gu>(x)  -  Qn(x)  >  y). 
In  short,  o)  =  u®  ^  (with  X  >  A)  yields  a  solution 
(i.e.  a  dynamically  possible  steady  flow)  in  the  entire 
domain  D  even  though  the  class  of  admissible  func¬ 
tions  is  restricted  to  those  supported  in  the  truncated 
domain  Da. 

We  do  not  have  an  existence  theory  in  D  without 
the  restriction  that  X  be  large  because  the  necessary 
a  priori  estimates  on  the  support  of  truncated 
solutions  are  available  (to  us)  only  for  large  X. 

All  of  the  statements  given  in  §4  concerning  the 
asymptotic  behavior  of  solutions  as  X  +  ®  carry  over 
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to  the  solutions  found  in  Theorem  5.2.  In  particular, 
the  location  Xq  of  the  limiting  point  vortex  is 
determined  by  (4.2)  and  hence  satisfies  the  conclusions 
of  Lemma  5.1.  That  is,  it  represents  an  equilibrium 
point  for  the  Hamiltonian  system 

||  -  -J7Hq(X)  . 

where  X(t)  describes  the  motion  of  a  (unit)  point 
vortex  moving  in  the  fluid  domain  D  in  the  presence 
of  an  irrotational  flow  with  streamf unction  -Qn.  In 
the  special  case  when  D  *  (i.e.  flow  past  a 

cylinder)  the  equilibrium  point  x  ■  x£  has  been 
computed  explicitly  by  FOppl  using  the  method  of 
images;  it  is  found  to  satisfy  the  equations 

2|xlx2  -  l xl 2  -  1,  1  -  2Qx2(1  -  lx|"4)  . 

§6.  VORTEX  PAIRS 

The  most  well-known  example  of  steady  vortex  flow 
in  R2  is  that  due  to  a  uniformly  translating  vortex 
pair.  Taking  account  of  the  odd  symmetry 
(wtXj^-Xj)  *  wfxj^Xj))  of  w  on  R2  this  problem 
corresponds  to  the  special  case  when 
D  *  Dq  =»  {x  e  R2  s  x2  >  0)  in  the  present  context. 
Furthermore,  solutions  also  always  possess  the  natural 
symmetry  condition 

(6.1)  wf-x^Xj)  *  u(xlfx2)  ; 

throughout  this  section  we  append  this  additional 
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requirement  to  the  definition  of  the  class  of 
admissible  functions  KX(DQ).  The  (energy)  functional 
Eq  is  defined  on  KX(DQ)  by 

(6.2)  E  («)  "  \  J  /  g0(x,  x*  )w(x)u»(x' )dxdx' 

D0  D0 

-  Q  /  x2w(x)dx 


where  the  Green  function  is 


(6.3) 


1  ^  X1 

90(x'x' }  =  47  109  - - 

<x. 


x|)2  +  (x2  +  xp2 
x{)2  +  (x2  -  x')2 


We  shall  exploit  the  especially  simple  form  of  the 

problem  on  Dg  in  order  to  give  sharper  results  than 

can  be  given  for  the  general  problems  of  §5. 

For  any  w  e  KX(DQ)  we  define  its  symmetrization  in 
* 

the  x1  variable,  u  e  Kx (Dg ) ,  in  the  usual  way;  namely, 
0)  (x1#x2)  »  <*»  (-xlfx2)  =  *u  (xx;x2)  (xx  >  0) 

where  *u(t;x2^  (a  monotonic  function  of  t)  is 

#u(t;x2)  *  meas.{x|  e  [0,®)  :  u>(xx,x2)  <  t)  . 

* 

We  observe  the  important  fact  that  Eq(o>  )  >  eq (“)f 
this  follows  since  g0(x,x*)  is  a  decreasing  function 
of  (Xjl  “  xj)  for  fixed  X2»x£  and  hence 


*  ★ 


/  /  gQ ( x, x ' )u>  (x)w  (x ' Jdx^dxj^ 


>  /  /  g0(x,x'  )u>(x)w(x'  )dxxdxx  , 


—00 
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a  consequence  of  standard  theory.  Therefore,  any 

maximizer  m  for  Eq  over  K^(DQ)  can  be  assumed  to 

* 

be  symmetrized,  u>  =  u»  . 

We  begin  our  analysis  of  vortex  pairs  by  stating 
an  existence  theorem  analogous  to  that  given  in  §5, 
providing  solutions  in  DQ  for  large  X. 

Theorem  6.1.  Let  a  =  (a1#a2)  be  fixed  so  that 
al,a2  >  2  1/4wQ} .  Then  any  (absolute)  maximizer 

“  “  “q  X  ^or  eQ  over  K*(D)  satisfies 
(6.4)  supp  «  c  (Ixjl  <  aj_  -  1,  0  <  x2  <  a2  -  1} 

provided  X  >  A 

where  A  is  a  sufficiently  large  constant  depending  on 
a,Q. 


Proof.  We  need  only  modify  the  proof  of  Theorem  5.2. 
By  virtue  of  condition  (6.1)  and  the  fact  that 

H0(X)  =  109  ji-  , 

the  only  critical  point  x  =  Xq  of  H0(x)  +  Qx2  is 
Xq  *  (0,l/4wQ).  As  before  we  find  that 
supp  m  c_  B^(Xq)  if  X  >  A  for  a  sufficiently  large 
constant  A,  and  this  yields  (6.4). 

The  above  theorem  ensures  that  the  maximizers 
u  =  UQfx  X  >  A  represent  solutions  of  the 

vortex  pair  problem  in  DQ.  These  solutions,  being 
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symmetrized,  then  have  the  form 
(6.5)  m  •  Xlfl,  8  ■  {x  6  Dq  *  G0w(x)  -  Qx2  >  y) 
with  fl  *  (Ixj^l  <  z(x2)>  for  a  function  z  defined 
implicitly  by  4»(z(x2),x2)  *  0.  It  is  straightforward 

it 

(using  «  -  w  )  to  verify  that  ♦  (x)  <  0  for 

i  2 

x^  >  0,  and  thus  z  e  CA  on  those  open  intervals 
where  z  >  0. 

The  requirement  in  Theorem  6.1  that  X  be  large 
is  more  restrictive  than  the  variational  theory  should 
require  and  is  imposed  mainly  for  technical  reasons.  On 
the  other  hand,  the  following  theorem  shows  that 
solutions  in  K^(DQ)  for  arbitrarily  small  X  are  not 
possible. 

Theorem  6.2.  There  does  not  exist  u>  *  u  e  K^(Dq) 
such  that  (6.5)  holds  (for  some  constant  y)  if  X  is 
sufficiently  small  depending  on  Q,  0  <  Q  <  •. 

In  the  proof  of  this  nonexistence  result  we 
require  three  lemmas  concerning  a  priori  estimation  of 
the  support  of  solutions  having  the  form  (6.5);  these 
lemmas  are  presented  first.  We  remark  that  in  (6.5)  it 
is  necessary  that  y  >  0;  otherwise,  the  containment 
(0  <  x2  <  |y|/Q)  c  n  contradicts  the  constraint 
meas.  n  =  1/X  <  «». 


Lemma  6.3.  For  all  x  e  supp  «  there  holds 
(6.6)  Q2|x1lx2  <  E(w)  »  i  /  b>(x)GQu(x)dx  . 


Proof.  Let  x  e  supp  o»;  we  observe  that 

(x£,x2)  e  supp  u  whenever  |x{|  <  IxjJ  by  virtue  of 

* 

the  symmetrization  u>  ■  u  .  Since  u  >  0  we  then  have 

x2 

Qx2  <  G0<o(x',x2)  <  f  |VG0«(x^,x^)|dx^  . 

Upon  integration  in  x{  this  yields 

2Q|x1lx2  <  /  | Vgqw( x * ) | dx  * 

I | < I x1 1 

o<x2  <*2 


<  {2\x.\x0fe  (J  |VGn»(x,)|2dx')1^ 
D0 

»  2{  |  Xj^l  XjEfw)}3^  , 


as  required. 


In  the  results  to  follow  we  make  use  of  the 
elementary  estimates 

(6.7)  g0(x, x* )  <  log(  1  +  p) 


C0  log  I'x-x1 1  if  >x~x’l  *  x2/2 


if  Ix-x'l  >  x,/2 


for  a  certain  absolute  constant  Cq. 
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Lemma  6.4.  For  all  x  6  supp  w  there  holds 
(6.8)  x2  <  CQ”1  max{l,log  fa  x2>  . 

Proof.  Let  Xwp^  ™  1  d®flne  P\  >  0*  If  *  x2 
then 


x2 

Gnw(x)  <  Cn  /  log  .  rrr  Xdx’  <  C  log(x2/p.) 

0  0  | x~x ' I <P x  ' 


If  Px  >  x2/2,  then 


X2 

G  w(x)  <  cn  /  log  r—  1 1  Xdx' 

0  0  |x-x'|<x2/2  |x  x  1 


+  Cnx,  /  |x-x'|“1Xdx' 

4  x2/2<|x-x* I <px 


<  CXx2px  . 

In  either  case  we  conclude 

Qx2  <  G0u(x)  <  C  max{l, log(x2/px)>  , 
and  so  (6.8)  is  proved. 


Lemma  6.5.  For  all  x  e  supp  w  there  holds 

q2  |  x  | 

(6.9)  |  Xi  I  <  CQ_1[l  +  ~l*!t  max {l, log  - — |]  . 

°  ^EU) 

Proof.  Recalling  (6.6)  we  define 
2 

y2  =  2E(u»)/q  Ixj^l  >  2x2;  then  whenever  x*  e  supp  w 
there  must  hold  |x{|  <  lx-,  1/2.  In 


with  xA  >  y2 


light  of  this  we  have 

(6.10)  Qx-j  <  Gn«(x)  <  /  gn(x,  x*  )to(x' )dx ' 

0  0<*^<y2 


I x^ | < | x1 I /2 


gQ(x,x'  )u)(x*  )dx'  . 


We  proceed  to  estimate  these  integrals  separately. 
First  terms 


+  C0x2  /  | x-x' |“1w(x' )dx* 

U  |x-x'|>x2/2 

°<x^y2 


<  CAx2  +  CAx2y2max{l,log(l/A^2  y2)> 

<  CAx2y2maxtl,  logd/A^2  y2)  >  ; 

here  (in  the  second  member)  we  observe  that  for 

|x  -  x'l  >  X2/2  there  holds  |x  -  x'l-1  <  Six  -  x'l-1, 

and  thus,  using  u  e  K^(Dq),  we  need  only  estimate 

/  I x-x ' | -1Adx* 

|x][-x^  |  <2max{y2, 1/Ay2) 

0<x’<y2 
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Second  term: 


|x^|<|x1l/2 


gQ(x, x' )w(x' )dx' 


<  Cnx-  /  |x-x'  |“1a»(x'  )dx* 

0  2  IxJMxJ/2 

<  Cx2 I  X^ I . 

Returning  with  these  estimates  to  (6.10)  we  obtain 
Qx2  <  Cx2[Xy2max{l,  log(l/l^  y2 ))  +  IxjJ-1]  ? 
the  desired  result  (6.9)  now  follows  by  the  definition 
of  y2. 


Proof  of  Theorem  6.2.  We  assume  that  0  <  A  <  1. 

Lemma  6.4  implies  that  there  is  a  bound  b2  *  b2(Q) 
such  that  x2  <  b2  for  x  e  supp  ».  With  b2  >  1 
fixed  we  may  estimate  the  energy 

0  <  E(«)  <  sup  Gn«  «  C(Q) A  log  2/A  ; 

a  u 

to  see  this  we  apply  the  reasoning  used  in  the  proof  of 

Lemma  6.5  (with  b2  replacing  y2)  to  get 

Gnw(x)  -  /  gn (x, x ' )u ( x ' )dx ' 

0«x'<b2 

<  CAx2b2max{l#  log(l/A^  b2 ) ) 

<  CAb2log(2/A1/2  b2)  . 

Lemma  6.5  now  implies  that  there  is  a  bound  b^  *  b^Q) 
such  that  IxjJ  <  b^  for  x  6  supp  w.  Obviously, 
these  bounds  (independent  of  A)  for  the  support  of 
w  together  give 
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/  u(x)dx  <  2Xb.b9  <  1 
D0 

provided  X  is  sufficiently  small  (depending  on  Q), 
and  this  is  a  contradiction  since  w  e  K^(Dq). 

In  view  of  the  above  nonexistence  result  it  is 
desirable  to  widen  the  class  of  admissible  functions  in 
order  to  be  able  to  give  a  variational  theory  without 
the  restriction  that  X  be  large.  We  therefore 
consider  the  new  class 

K'  (D.)  *  {<•»  t  /  wdx  <1,  0  <  w  <  X  a.e.  in  Dft)  , 

where  we  also  continue  to  impose  (6.1).  The  existence 

* 

of  an  absolute  maximizer  w  *  w  for  Eq  over 
k^(dq)  can  be  proved  for  arbitrary  0  <  X  <  •  using 
the  methods  of  this  paper.  The  support  of  any  such  w 
can  be  bounded  (depending  on  Q)  using  the  proofs  of 
Theorem  6.1  and  Theorem  6.2  respectively  for  the  cases 
of  large  X  and  small  X.  Also,  as  we  expect, 

■  1  X  large 

/  wdx  ' 

Dq  <  1  X  small  , 

although  a  sharp  estimate  for  the  value  of  X 
separating  these  two  cases  is  not  clear.  We  omit  the 
details  here. 

We  close  our  discussion  of  vortex  pairs  with  a 
description  of  an  alternate  variational  approach  to  the 
problem.  Rather  than  prescribing  the  speed  Q  (in  the 
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functional  Eq)  we  prescribe  instead  the  vortex 
impulse  P{«);  then  Q  must  be  determined  along  with 
the  extremal  u  as  the  Lagrange  multiplier  correspond¬ 
ing  to  the  vortex  impulse  constraint.  The  variational 
problem  is  then  to  maximize  the  energy  functional 

E(w)  -  i  /  «(x)Gn«(x)dx 

Do 

subject  to  the  constraint 

/  x0u(x)dx  ■  P  (0  <  P  <  •) 

Do 

over  the  class  of  admissible  functions  «  8  K^(Dq). 

The  existence  of  a  maximizer  w  •  u  ■  uip  .  can 
be  proven  for  arbitrary  0<A<«,  0<P<*»,  and  any 
such  solution  can  be  shown  to  have  the  form  (6.5)  for 
some  Q  >  0,  u  >  0.  This  alternate  approach  has  some 
particular  advantages  when  used  in  the  analogous 
(axisymmetric)  problem  of  vortex  rings,  and  has  been 
worked  out  in  detail  in  Friedman  and  Turkington  [2]y 
the  variational  formulation  itself  is  an  adaptation  of 
that  proposed  by  Benjamin. 

The  speed  Q  is  determined  by  <■>  as  follows: 

,  (x,+xl) 

(6.11)  Q  /  «o(x)dx  *  4-  /  /  — = — ^  u(x)«(x')dxdx' 

D  2"  D0  D0  lx-n2 

(x'  *  (x^-x^))  . 

To  show  this  we  observe  that  since  w  = 
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integration  by  parts  yields 
0  ■  /  w*  dx 

D0  2  -o 


/ 

Dr 


“(G0tt)x. 


dx  - 


Q  / 

Dr 


wdx 


thus,  noting  an  obvious  cancellation, 

Q  /  wdx  -  /  /  gQ(x,  x*  )w(x)o)(x*  )dxdx* 

Dq  Dq  Dq  2 


=  /  /  -J—  hn(x,  x'  )w(x)u»(x'  )dxdx'  , 


Do  Do  2 


which  is  (6.11). 

The  asymptotic  analysis  as  X  ♦  •  remains  valid 
for  the  solution  in  the  alternate  approach.  Now  we 
find  the  convergence  (in  the  distributional  sense) 

«p  X(x)  ♦  6(x  -  X*)  as  X  ♦  - 
where  Xp  *  (0,P)  by  virtue  of  the  vortex  impulse 
constraint  (and  (6.1)).  Also,  as  a  consequence  of 
(6.11),  we  have 

Qp,X  *  Qp  °  1/ 4»P  as  X  ♦  •  . 

Finally,  the  limiting  form  of  solutions  as  discussed  in 
§4  remains  equally  valid  for  x  as  X  ♦  ■». 


APPENDIX:  Vortex  energy  and  impulse 

We  include  this  discussion  in  order  to  clarify  the 
physical  basis  of  the  concepts  of  energy  and  impulse  of 
an  ideal  fluid  flow  defined  in  §1.  These  concepts  are 
developed  in  the  literature  only  for  the  motion  of  a 
body  through  an  irrotational  flow  or  for  a  vortex 
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motion  in  an  unbounded  fluid  (without  interior 
boundaries).  The  corresponding  results  in  the  more 
general  context  needed  here  -  the  motion  of  a  body 
through  a  rotational  flow  -  do  not  seem  to  be  available 
in  the  standard  references;  see,  for  example,  Batchelor 
Cl]. 

We  consider  a  body  D  (compact,  connected, 
simply-connected,  smooth  boundary)  moving  with 
translational  velocity  U(t)  and  zero  angular  velocity 
through  an  ideal  fluid  at  rest  at  infinity;  the  fluid 
occupies  the  exterior  domain  D  *  r2\d.  The  (fluid) 
velocity  field,  u  ■  u(x,t),  relative  to  a  frame  fixed 
in  D  is  governed  by  the  equations 
(A.l)  V»u  “0  in  D,  v*u  *  0  on  3D  , 

(A. 2)  Ut  +  U*Vu  »  -Vp  -  U  in  D  ; 

v  denotes  the  unit  normal  on  3D  directed  exterior 
to  D.  The  (fluid)  velocity  field  relative  to  a 
corresponding  frame  fixed  at  infinity  is  then  given  by 

A  A 

U  =  U ( X , t )  ■  u (x, t )  +  U(t)  . 

We  shall  require  that 

(A. 3)  |u|  =  0(|x|"2),  | Vu|  »  0(|x|’3)  as  lx|  ♦  -  ; 

this  restriction  -  adequate  for  our  purposes  -  is  more 
than  is  necessary  for  a  general  discussion  of  these 
concepts.  The  vorticity,  u  =  u(x,t),  is  given  by 
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The  dynamical  equation  (A. 2)  can  be  written  in  either 
of  the  following,  well-known,  equivalent  forms  1 
(A. 4)  ufc  *  wJu  -  Vb,  b  ■  j  |u|*  +  p  , 

(A. 5)  +  u*7w  ■  0  ; 

recall  J(a^,a2)  *  (a2,-a1)  . 

We  now  define  two  functionals  of  any  velocity 

A 

field  u  satisfying  (A. 3)  and  the  condition  that  e 
has  bounded  support  in  D.  The  fluid  energy  is  defined 
to  be 

(A. 6)  E *  «  i  J  I u| 2dx  . 

*  D 


The  fluid  impulse  is  defined  to  be 
(A. 7)  P*  -  (pJ,P2)  with 

f  it  * 

P 1  =  /  x-utdx  +  /  x9u*Jvds 
D  z  3D  z 

« 

#  * 

P2  *  -  /  x.wdx  -  /  x.u»Jvds  . 

D  1  3D  1 

The  physical  significance  of  E*  is  entirely  standard; 

that  of  P*  is  not  however.  In  an  appropriate  sense 

P*  takes  the  place  of  the  fluid  momentum  which  cannot 

be  defined  directly  because  the  integral 

/  udx 
D 
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is  divergent  according  to  (A. 3).  Indeed,  if  the  lack 
of  convergence  is  simply  ignored  then  the  expressions 
(A. 7)  can  be  derived  formally  from  the  latter  integral 
by  integration  by  parts. 

In  order  to  clarify  the  meaning  of  (A. 7)  we 
provide  the  following  proposition  in  which  the  impulse 
is  related  to  the  reaction  force  of  the  fluid  on  the 
body  D  when  accelerated. 

A 

Proposition.  Let  u  satisfy  (A.l,  A. 2,  A. 3)  and  the 
condition  that  «  has  bounded  support  in  D.  Then 

* 

(A. 8)  F  —  -  /  pvds  -  . 

3D  at 

We  shall  derive  the  first  component  of  (A. 8)  only, 
as  the  other  is  similar;  namely,  we  shall  show 

j  a 

(A. 9)  -  /  pv  ds  *  -Sr-  {/  x.wdx  +  /  x,u»Jvds}  . 

3D  1  D  3D 

Using  the  definition  in  (A. 4)  we  can  write 

-  /  pv  ds  =  /  i  |u|2v.ds  -  /  bv  ds  , 

3D  3D  3D 

and  so  we  consider  these  two  integrals  separately. 

First  integrals 

/  ^  |u|2v  ds  =  /  i  |u|2v  ds  -  /  u»Uv.ds 

3D  1  3D  *  1  3D  1 

*  /  (i  |u|2)dx  -  /  U'Uv.ds  . 

D  9X1  2  3D  1 
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0  ■,  *  2  *  *  * 

Applying  the  identity  (j  |u|^)  =  u»u2  +  7»(u1u)  we 

then  have 

/  i  |u|2v  ds  »  /  wu-dx  +  /  [ ( v»U)u.  -  (u*U)v  ]ds 

3D  Z  1  D  3D 


*  /  a»u,dx  +  U~{f  u»dx  +  /  u*Jvds} 
D  1  L  D  3D 


-  /  <uu,dx  ; 
D  * 


here  we  use  the  £act  that  the  (total)  circulation 

it  a 

(A. 10)  r  *  /  wdx  +  /  u«Jvds  =  0 

D  3D 


since  integration  by  parts  is  justified  according  to 

(A.3).  Continuing  we  calculate 

/  wu,dx  *  /  Vx0*(wu)dx 
D  Z  D  * 


=  -  /  x,u»Vwdx 
D  * 


=  /  x  w  dx  , 
D  Z  Z 


by  virtue  of  equation  (A. 5).  Therefore, 

(A. 11)  /  \  |u|2v  ds  *  /  x0u»dx  . 

3D  dt  D  1 
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Second  integral) 


-  J 

3D 


bv^ds 


/ 

3D 


3x, 

4 

3  T 


ds 


r  3b 

'  iD  2  77 


/  x2(u  - 
3D  *  Z 


J  X,U  *TdS 
3D  *  * 


ds 


Uju) •tds 


with  t  ■  Jv  (unit  tangent)  on  3D;  here  we  use 
integration  by  parts  around  3D  followed  by  equation 
(A. 4).  Therefore, 

A  * 

(A. 12)  -  /  bv  ds  -  /  x,u»Jvds  . 

3D  at  3D  * 

The  required  identity  (A. 9)  now  follows  by  combining 
(A. 11)  and  (A. 12). 


The  terra  "impulse"  is  derived  from  another 

*  * 

interpretation  of  P  .  We  say  that  u  is  expressed  in 
terras  of  impulsive  pressures,  p*,  and  impulsive 
(body)  forces,  f*,  if  it  can  be  written 
(A. 13)  u  =  -Vp  +  f  with  |f*|  ■  0(|x|""2~e) 

as  I  x  |  ♦  »,  e  >  0  . 

It  then  follows  that 

(A. 14)  P*  =  -  /  p  vds  +  J  f  dx  , 

3D  D 
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and  this  expression  is  independent  of  the  choice  taken 
in  (A. 13).  A  further  discussion  is  given  in  Lamb  [3] 
§1.11. 

In  the  case  that  the  body  is  assumed  to  be 
symmetric  about  the  axis  and  the  flow  likewise  the 

expressions  (A. 6)  and  (A. 7)  take  an  especially  simple 
form;  this  is  the  situation  considered  in  §5.  In 
keeping  with  the  notation  established  earlier  (see  §5) 
we  shall  now  let  D  and  D  denote  respectively  the 
intersection  of  the  body  and  the  (fluid)  domain  with 
the  upper  half  plane  (x2  >0}.  We  assume  also  that 
U  a  q ( 1 , 0 ) ,  0  <  Q  <  «*.  These  symmetry  assumptions 
(along  with  (A.3))  imply  that  u>(x)  ■  -u(x), 

a  a  a  a  _ 

uL(x)  *  u1(x),  u2(x)  »  -u2(x)  for  x  -  (xj^-Xj)# 

x  e  D.  Thus,  we  find 

(A.  6 '  )  4  E *  »  4  J  I  ul  2dx  , 

z  z  D 

(A. 7')  i  P*  -  /  x,udx  +  /  x~u*Jvds,  *  0  • 

1  1  D  *  3D  ^ 

We  claim  that  these  quantities  can  be  expressed  in  the 
alternate  forms: 

(A. 15)  j  E*  -  /  wGwdx  +  j  m2  , 

(A. 16)  4  P?  *  /  nwdx  +  MQ,  P?  «  0  , 

2  1  D 

A# 

where  the  induced  mass  of  D  is  defined  by 


(A. 17) 


Ms/  | V ( n  —  X2 ) I 2dx  ; 

D 

we  retain  the  definitions  of  G  and  n  from  §5. 

In  the  derivations  of  these  identities  we  let 
n  *  n  -  X2?  then  the  streamfunction  for  the  velocity 

*  A  A 

field  u  is  given  by  ^  =  Gut  -  Qn. 

Derivation  of  (A. 15): 

4  E*  =  4  /  I vj | 2dx 

2  D 

-  4  /  |7Gu|2dx  -  Q  /  7Gw7ndx  +  4  Q2  /  |7^|2dx 

D  D  2  D 

-  4  /  “Gudx  +  4  Q2  /  I 7n| 2dx  . 

D  D 

Derivation  of  (A. 16): 

A 

4  P -i  *  /  nwdx  +  /  nA^dx  -  /  n  |-i  ds 
2  1  D  D  3D  d 

A 

=  /  nwdx  -  /  ♦  4^  ds 

D  3D 


■  /  nwdx  +  Q  f  n  4^  ds 
D  3D 


=  /  nwdx  +  Q  /  |Vnl2dx  . 

D  D 

As  a  consequence  of  these  identities  we  see  that 
the  energy  and  impulse  may  both  be  split  into  a 
"rotational  part"  and  an  "irrotational  part"?  the 
functionals  E(u>)  and  P(w)  defined  by  (1.10)  and 


-39- 


(1.11)  represent  the  rotational  parts  respectively, 

1  o 

while  MQ  and  MQ  with  definition  (A. 17)  represent 

the  irrotational  parts  respectively. 
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